Molecular Orbital Shift of Perylenetetracarboxylic-Dianhydride 
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Using low-temperature scanning tunneling microscopy we find that perylenetetracarboxylic- 
dianhydride on Au(788) exhibits three coexisting adsorption phases. Single-molecule tunneling 
spectroscopy reveals orbital energies, which differ in the different adsorption phases. Density func- 
tional theory calculations associate the experimentally observed submolecular corrugation to the 
spatial distribution of the second- to- lowest unoccupied molecular orbital. We tentatively attribute 
the orbital shifts to a varying number of hydrogen bonds. 

PACS numbers: 68.37.Ef,73.20.At, 73.61. Ph 



Interest in organic molecular films is fueled by techno- 
logical advances in organic devices 0, Q , potential appli- 
cations in molecular switches 0, Q , and functionalization 
of surfaces [1, 0, R as well as by the fundamental under- 
standing of the film-metal system. Typically, the geo- 
metric and electronic structure of the molecule-substrate 
interface is of substantial importance. Therefore consid- 
erable effort has been made to investigate this interface 
and to maximize its homogeneity 8] . Owing to the com- 
plexity of molecular shapes and interactions, however, 
perfectly homogeneous interfaces usually remain elusive 
and coexistence of structural domains has been reported 
for various molecular systems 9]. Detailed characteri- 
zation of the variability of the electronic structure on a 
single molecule scale requires high lateral resolution and 
has, consequently, hardly been achieved 

A model system of molecular electronic studies at 
the molecular level is 3, 4, 9, 10-perylenetetracarboxylic- 
dianhydride (PTCDA). The adsorption geometry of 
PTCDA on metal surfaces was analyzed previously as ex- 
cellent lyreviewed by Witte and Woll and Barlow and 
Raval [13 • Adsorption characteristics on Au(lll) de- 
serve special interest since this surface exhibits a 22 x 
reconstruction. Chizhov et al. found that room- 

temperature adsorption leads to two coexisting adsor- 
bate geometries, namely a herringbone structure and a 
square phase, where four PTCDA molecules adsorb on 
the boundary of a square with equivalent axes of adja- 
cent molecules including an angle of 90°. 

Here, we report on low- temperature scanning tunnel- 
ing microscopy (STM) of PTCDA on Au(788). We ob- 
serve three coexisting PTCDA superstructures on this 
surface. Single-molecule tunneling spectroscopy reveals a 
peak above the Fermi level whose energy strikingly varies 
between these geometric superstructures. According to 
density functional theory (DFT) calculations for the 
molecule-surface complex the experimentally observed 
spectral feature is predominantly due to the second-to- 




FIG. 1: Constant-current STM image of 0.25 ML PTCDA on 
Au(788) (310 A X 310 A, y = 0.75 V, / = 200 pA). Crystal- 
lographic directions along and perpendicular to step edges are 
indicated. A, B^, and C denote different PTCDA adsorption 
phases. Inset: clean surface area with v-shaped discommensu- 
ration lines. Left hand side: sketches of the observed PTCDA 
superstructures including lattice vectors. 



lowest unoccupied molecular orbital (LUMO+1). This 
interpretation is in some contrast to previous work for 
PTCDA on other metal surfaces where experimental fea- 
tures have been attributed to the LUMO on the basis of 
calculations for the free molecule 0, Q, . We argue 
that the observed shift of the orbital energy is most likely 
caused by hydrogen bonds. 

We used a home-built ultrahigh-vacuum (UHV) low- 
temperature STM with a base pressure of better than 
1 X IQ-^Pa. The Au(788) sample and W tips were 
cleaned by argon ion bombardment and annealing. Crys- 
talline order was checked by low-energy electron diffrac- 
tion, and surface cleanliness was monitored by STM. 




FIG. 2: Sequence of constant-current images recorded at characteristic sample voltages. The top (bottom) part of each image 
shows B-domain (A-domain) molecules (120 A x 120 A, / = 200 pA). A periodic grayscale is used to highlight the contrasts on 
each terrace. 



PTCDA was evaporated onto the Au(788) surface held 
at room temperature from a degassed tantalum cru- 
cible keeping the pressure of the UHV recipient below 
1 X 10~^Pa. The rate of deposition, monitored by a 
quartz microbalance, was 0.5 ML min~^. One monolayer 
(ML) of PTCDA is defined here as one molecule per 
100 A (corresponding to the density of the herringbone 
superstructure). The prepared sample was then cooled 
down to 4 K within 30 min. For detection of the differen- 
tial conductivity {dl/dV) we applied a standard lock-in 
technique. The tip state was checked by spectra of the 
electronic surface state of clean Au(788) acquired before 
and after measuring the dl/dV signal of single PTCDA 
molecules. 



Density-functional calculations were performed with 
the code Siesta \1& using a double-^ basis set plus po- 
larization functions [13 for the valence electrons and the 
Generalized Gradients Approximation (GGA) for the 
exchange-correlation functional. The Au(lll) surface 
serves as a model of the terraces of Au(788). Various 
supercell sizes in a slab geometry were explored to take 
into account different packing densities of PTCDA super- 
structures. For simulating experimental spectra of the 
differential conductivity we calculated the local density 
of states at the tip apex being 7 A above the surface. 
Due to the weak molecule-surface interaction local and 
semi- local approaches of DFT fail to yield the correct ad- 
sorption parameters of PTCDA on Au(lll). Here, we are 
interested in a qualitative analysis of the electronic struc- 
ture. A recent incidence x-ray standing wave experiment 
performed on Ag(lll)-PTCDA determines the vertical 
distance of the perylene core to 2.8A Hence, we 
have frozen the height of the molecular atoms at 2.8 A 
and allowed a full relaxation of the other molecular co- 
ordinates in order to find the correct packing geometry 
for the two molecules constituing our supercell. Our cal- 



culations show that increasing the adsorption height by 
1 A leads to a basically free molecule, which slightly in- 
teracts with the surface. Decreasing the molecule-surface 
distance by 0.3 A leads to a very large mixing of molec- 
ular orbitals because the molecule is already probing the 
strongly repulsive substrate potential. 

Figure Q shows a constant-current STM image of 
0.25 ML PTCDA on Au(788). Crystal directions along 
ascending steps, [211], and parallel to the step edges, 
[Oil], are indicated. The inset shows a small-scale STM 
image of a clean area to highlight the discommensura- 
tion lines. These lines of higher corrugation are stacking 
fault domains separating face-centered cubic (fee) and 
hexagonal close packed (hep) surface layer stacking ar- 
eas, similar to the case of the 22 x ^/3 reconstruction on 
Au(lll). Au(788) is thus patterned in two dimensions: 
terraces with a 38 A periodicity along [211] and discom- 
mensuration lines with a periodicity of 72 A along [OlT|. 
We find that the Au(788) reconstruction is not perturbed 
by PTCDA adsorption. A similar observation has been 
reported previously for the 22 x >/3 reconstruction on 
Au(lll) [1^. 

Three adsorption phases of PTCDA coexist on 
Au(788) (labelled A, Bn, and C in Fig. [ij. In domain 
A we find the herringbone structure with lattice vector 
dimensions a = (11.2 ± 0.7) A and b = (17.9 ± 0.9) A 
(sketches on the left hand side of Fig. Three preva- 
lent orientations of the shorter unit cell edge with respect 
to [211] are observed, namely (33 ±4)°, (50 ±4)°, and 
(77 ±3)°. For Au(lll), a different set of orientations 
(22°, 39° and 55°) was reported j20| and rationalized in 
terms of a minimum misfit between the molecule film and 
the substrate. On a stepped surface like Au(788) this ef- 
fect appears to be suppressed by the possibility to align 
the molecules with the step edges. 

Superstructure B is a square phase with a lattice con- 
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FIG. 3: Experimental (crosses) and calculated (full line, fitted 
to the experimental maximum, the peak includes a numerical 
broadening of 0.25 eV) dl/dV spectra of a single PTCDA 
molecule in the herringbone adsorption domain (A) . Feedback 
loop parameters: / = 200 pA, V = 1.4 V. 

stant of a = b = (14.4 ± 1.2) A. Equivalent axes of ad- 
jacent molecules are rotated by 90° with respect to each 
other. Several rotational domains, labelled Bi, B2, and 
B3 in Fig. n are observed. Chizhov et al. reported a 
similar superstructure on Au(lll) fl^. 

Adsorption phase C differs from structure B by the end 
of a molecule pointing towards the center of an adjacent 
molecule. 

We now turn to the electronic properties of the 
PTCDA molecules in the superstructures. Figure [2 
presents STM images recorded at three characteristic 
voltages. The images show four PTCDA covered Au(788) 
terraces where A and B domains are visible in the lower 
and upper parts, respectively. In B domains, circular pro- 
trusions are frequently observed between four PTCDA 
molecules. We attribute these protrusions to single 
trapped Au atoms in analogy to single Fe atoms sur- 
rounded by trimesic acid molecules l2l|. At 400 mV all 
PTCDA molecules are imaged as structureless rectan- 
gles. For 800 mV topographic changes of the A phase 
are hardly visible while intramolecular structure emerges 
in the B phase, the former rectangular shape of the 
molecules being lost in favor of a hexagonal shape. Fi- 
nally, at 1200 mV submolecular structure is observed for 
both domains. These observations strikingly demon- 
strate that the electronic structure underlying the in- 
tramolecular contrast depends on the adsorption phase 
of PTCDA molecules. 

To interpret the data we first trace the origin of 
submolecular structure back to the involved molecular 
orbit als. Figure |31 shows the differential conductivity 
dl/dV as acquired in the center of a single PTCDA 
molecule embedded in the herringbone adsorption do- 
main. We find a single peak centered at ~ 1.1 eV. The 
calculated dl /dV spectrum, shown as a line, matches the 
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FIG. 4: Spectra of dl /dV vs. sample voltage on bare Au(788) 
(dots), on B-domain PTCDA (open circles), and on A-domain 
PTCDA (solid circles). Maxima of dl/dV occur at ^ 0.8 V 
and ^ 1.1 V for B- and A-domain PTCDA, respectively. The 
feedback loop was frozen at / = 200 pA, V = 850 mV. Data 
for B-domain PTCDA was multiplied by 4. 



experimental spectrum. In the calculation, the main con- 
tribution to the conductivity originates from the density 
of states of the LUMO+1. The LUMO+2 is basically de- 
generated with the LUMO+1 and also contributes to the 
spectroscopic signature giving a certain broadening of the 
peak. Since the intramolecular structure for herringbone- 
phase molecules appears above 1 V we ascribe most of 
the submolecular structure to the spatial distribution of 
the LUMO+1. 

We now return to the differences between domains A, 
B, and C observed in Fig. [21 Intramolecular structure for 
the square adsorption domain B appears at lower volt- 
ages than for the herringbone phase. Molecules at a do- 
main boundary surrounded from either side by molecules 
in different arrangements and likewise molecules in ad- 
sorption phase C reveal submolecular structure for in- 
termediate voltages. Obviously, the energetic position of 
the LUMO+1 depends on the local environment of the 
individual molecule and thus differs for molecules from 
different domains. We further tested this hypothesis 
by single-molecule spectroscopy of the differential con- 
ductivity acquired above the center of B and A phase 
molecules (Fig. 0| open and solid circles, respectively). 
The spectral signature of B phase PTCDA is located at 
~ 0.8 eV, i.e. shifted by ~ 0.35 eV compared to the 
A phase. Owing to the similar intramolecular contrast 
in both phases we can safely attribute this peak to the 
LUMO+1. Since the feedback loop was frozen at 850 mV 
the LUMO+1 contribution of the A-phase molecule has 
not yet developed leading to a smaller tip-molecule dis- 
tance than in the case of the B-phase molecule. To show 
both peaks on the same scale we multiplied the data ac- 
quired on the B-phase molecule by a factor of 4. Spectra 
of occupied electronic states, i. e., at negative voltages 
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(—2.5 V <V<0) turned out to be featureless. Our cal- 
culations predict, for instance, the LUMO of PTCDA to 
be pinned at the Fermi level 19]. The spatial distribu- 
tion of the LUMO density of states is found to adhere 
strongly to the molecule and to be less evanescent into 
vacuum than the LUMO+1. This leads to its efficient 
suppression from the dl/dV spectra (as observed in our 
experiments). 

The calculations indicate that apart from the molecule- 
substrate interaction the inter molecular one plays a role. 
The latter is mediated by hydrogen bonding between ad- 
jacent molecules which is poorly accounted for by our 
numerical approach. Nevertheless, our calculation for 
free standing molecules at the adsorbed geometries yields 
a widening gap between molecular levels. As a conse- 
quence the packing density and the actual adsorption 
geometry influence the energetic position of molecular 
orbitals. The ability of DFT to reproduce workfunction 
changes even in the case of physisorbed molecules j22j . 
led us to compute and compare the workfunction changes 
with the measured ones. The experimental change is of 
the order of 0.1 eV which agrees with the change of the 
two molecular phases frozen at 2.8 A from the surface 
atomic plane, implying that the molecules lie at the same 
molecule-surface distance for both phases. However the 
d//dV calculations show small differences between both 
phases and are unable to explain the significant experi- 
mental orbital shift. Different combinations of molecular 
heights for the two phases led to inconsistent results and 
to disagreement with the measured workfunction change. 
Hence, we conclude that the orbital shift is not due to 
a different molecule-surface distance for the two phases. 
Other effects such as the influence of the tip's electric 
field have been included in the calculation by computing 
the electronic structure for different values of external 
electrical fields without a sizeable effect on the dl/dV 
peak position. Finally, the modification of the metallic 
surface is unlikely due to the very weak physisorption in- 
teraction. These facts lead us to favor H-bond mediated 
molecule-molecule interaction, where the role of the sur- 
face may not be negligible due to the important change 
of charge transfer that the orbital shift implies. 

In summary, we found that PTCDA exhibits three 
distinct coexisting adsorption phases on Au(788). The 
adsorption domains can be distinguished by their geo- 
metric arrangement and by their electronic properties. 
By performing DFT calculations the experimentally sub- 
molecular corrugation is ascribed to the spatial distribu- 
tion of the LUMO+1 density of states. The shift of the 
LUMO+1 energy as seen in single-molecule spectroscopy 
measurements in the distinct superstructures might be 
ascribed to the different PTCDA packing densities and 
to the particular adsorption geometries. 
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